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Abstract. Exciton condensation in a two-band Hubbard model on a square lattice is studied 
with variational Monte Carlo method. We show that the phase transition from an excitonic 
insulator to a band insulator is induced by increasing the interband Coulomb interaction. To 
examine the character of the exciton condensation, the exciton pair amplitudes both in fc-space 
and in real space are calculated. Using these quantities, we discuss the BCS-BEC crossover 
within the excitonic insulator phase. 


1. Introduction 

Exciton condensation is one of the most interesting phenomena in strongly correlated electron 
systems and it has been extensively studied since 1960s mmm- The exciton is a bound electron- 
hole pair mediated by the (attractive) Coulomb interaction between them. When the binding 
energy of an exciton exceeds the energy gap between the electron and hole bands, the system has 
an instability towards condensation of excitons. Despite the extensive studies for half a century, 
the example of exciton condensation which is generally accepted is quite limited. One of the 
reasons is that the exciton is a neutral quasiparticle and difficult to be detected in experiments. 
Therefore, any conclusive evidence for exciton condensation in real materials is essential for 
further progress. 

Among the limited candidates, the transition metal dichalcogenide of TiSe 2 is one of the 
most notable materials for exciton condensation. It shows a commensurate CDW transition at 
T c ~200K with small structural distortion but the Fermi surface partially remains due to the 
imperfect opening of the energy gap. The origin of the CDW transition is still controversial 
and one of the possible mechanisms is exciton condensation mu- The observed flat band 
structure just below the Fermi energy [6] is characteristic of exciton condensation, although it 
is not a direct evidence. Another possible mechanism for the CDW transition is a Jahn-Teller 
effect which results from the electron-phonon interaction [7j. Recently, it is proposed that both 
mechanisms work cooperatively for the CDW state mm- In addition to the CDW phase, TiSe 2 
exhibits superconductivity with dome-shaped transition temperature induced by intercalating 
Cu atoms between layers [TOj or by applying pressure m ■ The origin of the superconductivity 
has not been clarified yet and the relation between the CDW and the superconductivity has 
attracted much interests in the context of the quantum critical point. 


Most of the theoretical studies for exciton condensation are limited to an ideal case, e.g., 
considering models in one dimensional, ignoring an intraband Coulomb interaction, a long- 
range Coulomb interaction, and a spin degrees of freedom. Although these models are useful to 
describe exciton condensation alone, it is generally difficult to discuss the competition between 
other competing phases such as magnetism and superconductivity due to the lack of realistic 
condition. Our aim is to establish a calculation scheme which can describe such a competition 
and to clarify the emerging mechanism of exciton condensation and superconductivity in low 
carrier density systems such as TiSe 2 - The variational Monte Carlo (VMC) method is one of 
the powerful methods for a strongly correlated electron systems and is suitable to this kind of 
problem. In this paper, we study a two-band Hubbard model with intra- and interband Coulomb 
interactions to investigate the detailed character of exciton condensation. The VMC method 
is used for the calculation of physical quantities in the ground state. We show that the phase 
transition from an excitonic insulator to a band insulator is induced by increasing the interband 
interaction. The importance of the intraband interaction is also discussed. Furthermore, the 
BCS-BEC crossover within the excitonic insulator is studied by calculating the exciton pair 
amplitudes both in fc-space and in real space. 


2. Model and method 

We consider a two-band Hubbard model on a two-dimensional square lattice defined by the 
following Hamiltonian, 
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where at (a = a, b) denotes the creation operator of an electron with wave vector k and 
spin a (=1\ i) on band a. U aa is an intraband Coulomb interaction and U ab is an interband 
Coulomb interaction between bands a and b. nf a = a\ a oiia denotes the number operator and 
nf = Y^(r n ?a- We set the non-interacting band energy as £ a (k) = 2t a (cosk x + cos k y ) and 
£b{k) = 2tb(cosk x + cos k y ) + Eq. The Fermi surface and the energy dispersion are shown in 
Figs. DJa) and (b), respectively. The electron and hole Fermi surfaces are perfectly nested with 
wave vector Q = (n, n) and have instability towards exciton condensation. 

Here, we introduce the following Gutzwillei'-Jastrow type trial wave function: 

|T> = P 3c Pg ] |$) . (2) 

The one-body part 1$) is obtained by diagonalizing the mean-field Hamiltonian, 
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Since we fix t a = t as an energy unit, e a {k ) is unchanged through the VMC calculation. On the 
other hand, tb and Eq in e&(fc) are optimized so as to minimize the variational energy. A k is a 
variational parameter for the exciton condensation and generates the spontaneous hybridization 
between bands a and b with wave vector Q = (ir, n). Although there are several types of exciton 
condensation depending on charge and spin degrees of freedom, they are energetically degenerate 
for the Hamiltonian without spin-dependent interactions such as Hund’s coupling. Therefore, 
we assume a spin-independent variational parameter as follows, 
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Figure 1. (color online) (a) Fermi surfaces and (b) energy dispersions of the noninteracting 
tight-binding energy band with electron density n = 2 . (c) Ground state phase diagram in the 
U aa /t — U ab /t plane. PM, El, and BI denote paramagnetic metal, excitonic insulator, and band 
insulator, respectively. PM is stabilized only for U ab /t = 0 due to the perfect nesting condition. 


A is an amplitude of the variational parameter and A controls the extent of exciton in fc-space. 
We have found that the exciton condensation can be described without A (i.e., A& = A) but the 
introduction of A greatly improves the variational energy and gives a better trial wave function. 
The Gutzwiller operator 
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in |T) is the one extended for the two-band system, i is a site index and 7 represents possible 
electron configurations at each site, namely, |0) = 10 0), |1) = |0t), 115) = Itlti)- 

The variational parameters gAs vary from 0 to 1, which control the weight of each electron 
configuration. Here, we classify the possible 16 local electron configurations into 10 groups by 
the local energy and set the same value of s for electron configurations with the same energy. 
The explicit grouping is shown elsewhere. 

The remaining operator 
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in |T) is the charge Jastrow factor, which controls the long-range charge correlations. Since 
the on-site Coulomb interactions U aa , U bb , and U ab affect not only the on-site but also the 
off-site (long-range) charge correlation, the charge Jastrow factor is essential for describing the 
appropriate trial wave function m- The band dependence is fully considered, namely, vf?, v bb , 

and vf b are independently optimized. For the spatial dependence, we assume that depends 
only on the distances (not on the direction), = v a ^(\r{ — rj |), and consider the range of 
R = |r* — rj | < L/ 2 for a square lattice of N = L x L. 

These variational parameters, A&, are simultaneously optimized so as to minimize 

the variational energy. For the optimization, we employ the stochastic reconfiguration 
method [T3J which works quite effectively for the VMC method. 


3. Results 

First, we study the stability of the exciton condensation and its competing phase by changing 
the intra- and interband interactions. The system size is 28x28 through the paper. Here, we 










Figure 2. (color online) Electron-hole pair amplitude in fc-space 4>(k) for (a) U ab /t = 4.6 and 
(b) U ab /t = 5.2. kp denotes the original (noninteracting) Fermi wave vector shown in Fig. UK a) 
with a red circle, (c) U ab /t dependence of the effective energy gap Eo,/t. Dotted line corresponds 
to the original value Ec/t = 6 shown in Fig. [[Kb). In all figures, U aa /t = U bb /t = 8.0. 


set the intraband interactions U aa = U bb for simplicity. Figure HKc) shows the obtained phase 
diagram in the U aa jt — U ab jt plane. The order parameter of exciton condensation is defined 
as A = a (b\ + Q a a ka + H.c.) and the finite A leads to a broken translational symmetry 
with ordering vector Q = ( n,ir ). (Note that A is different from the variational parameter A& 
defined in eq. (4).) Since the Fermi surfaces of bands a and b are perfectly nested as shown in 
Fig.[D(a), an infinitesimally small U ab induces the exciton condensation and the system becomes 
fully gapped. This phase is called excitonic insulator which has been extensively studied so 
far. As for the trial wave function, we have confirmed that the charge Jastrow factor Pj c is 
quite important for the stability of exiton condensation through lowering the variational energy 
of excitonic insulator phase. This is because the exciton pair has nonnegligible long-range 
component (at least several lattice constant) even with on-site Coulomb interaction alone. The 
details are discussed later in Fig. [21 When U ab increases further, the bands a and b are completely 
decoupled and the band insulator is realized. In the band insulator phase, the band a is fully 
filled ( n a = 2) and the band b is empty ( rib = 0) to avoid the energy loss of U ab despite the energy 
loss of U aa . Therefore, the region of the excitonic insulator becomes larger as U aa increases m- 
This phase diagram is consistent with the previous variational cluster approximation (VCA) 
study both qualitatively and quantitatively m • It is interesting to note that the excitonic 
insulator phase is restricted to a quite narrow region when the intraband interaction is absent 
( U aa = 0). This indicates the importance of the intraband interaction to stabilize the exciton 
condensation. Since most of the previous studies for exciton condensation neglect the intraband 
Coulomb interactions, its stability might be underestimated. 

Next, we calculate the electron-hole pair amplitude in fc-space defined as, 

H k ) = ^2( b k+Qa a kcr + H.c.) (7) 
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Figure [2] shows the behavior of 4>(k) in the first Brillouin zone. For U ab jt = 4.6, 4>(k) has a sharp 
peak at k = kp, which is an effective (or renormalized) Fermi wave vector in the wave function 
14/). A sharp peak structure in fc-space corresponds to a widely extended electron-hole pair in 
real space. The dotted line in Fig. EJa) corresponds to the original (noninteracting) Fermi wave 
vector kp and we can see that \kp\ > |fcp|. This is because the effective energy gap Eq in |4/) 
is reduced from the original value ( Eo,/t < Eq/1 = 6) and the Fermi surface is enlarged. On 
the other hand, as shown in Fig. [2Kb), the peak structure is rather broad for U ab ft = 5.2 and 
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Figure 3. (color online) Electron-hole pair amplitude in real space <f>(R) and R 2 \cj)(R)\ 2 for (a) 
U ab /t = 4.6 and (b) U ab /t = 5.2 at U aa ft = U bb /t = 8.0 (same with Fig. E]). 


it indicates a tightly bounded electron-hole pair in real space. The effective Fermi wave vector 
is smaller than the original one, |fcp| < |fcp|, and the effective energy gap is larger than the 
original one, Eg ft > Eg ft. The U ab dependence of Eg ft is shown in Fig. E]c) and this behavior 
is explained as follows. The character of the exciton condensation is determined by the balance 
between intra-(17 aa ) and interband ( U ab ) interactions. When the effect of U aa is dominant 
( U ab /t < 4.8), the effective energy gap is renormalized to a smaller value, Eg ft < Eg ft , to 
decrease n a and avoid the loss of U aa . As U ab increases for fixed U aa , Eg ft increases and 
exceeds the original value around U ab /t 4.8. For U ab /t > 4.8, the effect of U ab is dominant 
and the effective energy gap is renormalized to a larger value, Es/t > Es/t, to decrease rib 
and avoid the loss of U ab . Finally, the bands a and b are completely decoupled and the band 
insulator is realized for U ab /t > 5.24. Note that we cannot observe the limit of the Bose- 
Einstein condensation (BEC) of excitons where the effective Fermi wave vector kp vanishes, 
namely, Eg ft > 8 (no overlap between the bands a and b ) and A / 0. However, the behavior 
of the variational parameters drastically changes across U ab /t ~ 4.8 and we expect that there 
is a BCS-BEC crossover around this U ab /t. In the VCA study, the BCS-BEC crossover occurs 
around U ab /t ~ 4.0 and the vanishing of kp is observed at U ab /t = 5.2 [15]. Although the VMC 
study depends on the form of trial wave functions used, it can treat large enough system sizes, 
as large as N = 28 x 28, to describe the bound electron-hole pair of exciton at least in the BEC 
region. The electron correlation is exactly treated within a 2x2 cluster in the VCA, but our 
results suggest that the cluster size dependence should be carefully considered. The difference 
between the VMC and VCA results should be clarified to further understand the BCS-BEC 
crossover and it is precisely our future problem. 

Next, to examine the character of exciton in real space, we calculate the electron-hole pair 
amplitude 4>ij defined as 

faj = ^2( b lr a jo) + H.c. (8) 
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Although in general 4>ij depends not only on the distance between r* and Tj, R = |r* — rj |, but 
also on the direction between r* and rj, we have found that the direction dependence of (j)ij is 
small at least for parameters studied here. Therefore, we assume 4> t j = (p(R) in the following. 
Figure 0 shows the behavior of <f>(R) for U aa /t = U bb /t = 8.0. For U ab /t = 4.6 (expected to 
be BCS type), (p(R) gradually decreases with increasing R and exhibits oscillatory behavior. A 
similar oscillatory behavior can be observed in the electron-electron pair amplitude ( c !t c ]t)> of the 
BCS superconductivity, which is proportional to sin(kpR) //?,. This is because the wave function 










of exciton condensation is composed of a coherent electron-hole pair and essentially identical 
to that of the BCS superconductivity with a coherent electron-electron pair [2], On the other 
hand, for U ab /t = 5.2 (expected to be BEC type), in the vicinity of the phase boundary between 
the excitonic and band insulators, 4>(R) decreases more rapidly and the oscillatory behavior is 
hardly observed. As U ab increases, the attraction between electron and hole increases and the 
radius of exciton becomes smaller in real space. These behaviors are consistent with the fc-space 
picture: a widely-extended (tightly-bounded) exciton in real space corresponds to a localized 
(an extended) exciton in fc-space. 

Finaly, we calculate the coherence length of electron-hole pair defined as 
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The coherence length is often used to distinguish the BCS region (£ > 1, weak coupling) from 
the BEC region (£ < 1, strong coupling). In Fig. [3l we plot the quantity R 2 \4>(R)\ 2 appearing in 
Eq. ((HI) . For U ab /t = 5.2, this quantity is almost negligible for R > 4, indicating that the system 
size (IV = 28 x 28) is large enough to describe the bound electron-hople pair of exciton. The 
estimated coherence length £ is 0.77. On the other hand, for U ab /t = 4.6, i? 2 |</>(i?)| 2 gradually 
decreases with increasing R, showing oscillatory behavior, and has significant contribution 
for rather large R. We have also found that the size dependence of the calculation becomes 
severe around this U ab /t and below. This indicates that there is a BCS-BEC crossover around 
this U ab /t, below which the coherence length rapidly increases in the BCS region where even 
IV = 28 x 28 is not large enough to estimate £. Although it is difficult to discuss the BCS-BEC 
crossover with and £ shown here, we can visualize the spatial distribution of exciton using 
these quantities, which helps for further understanding the exciton condensation. The accurate 
estimation of £ and other physical quantities in real space is an important future problem. 

4. Summary and discussion 

In this paper, we have studied the two-band Hubbard model with perfectly nested electron 
and hole Fermi surfaces to discuss the stability of the exciton condensation with changing the 
Coulomb interactions. We have shown that the phase transition from an excitonic insulator to 
a band insulator occurs by increasing U ab and the region of the excitonic insulator is extended 
for larger U aa . We have also calculated the exciton pair amplitude both in fc-space and in real 
space and proposed that the BCS-BEC crossover occurs around U ab /t ~ 4.8 for U aa /t = 8.0. 

The extension of the VMC study for more realistic models is straightforward. We have 
checked that the introduction of a long-range Coulomb interaction greatly extends the region 
of the excitonic insulator phase in the phase diagram [16]. This is because the long-range 
Coulomb interaction enhances the long-range exciton pairs and leads to a more stable excitonic 
insulator phase. In such a case, the charge Jastrow factor Pj c becomes more important 
because the screening effect of long-range Coulomb interaction, which is an essential element 
for the exciton condensation, will be properly included through Pj c . On the other hand, the 
exciton condensation is greatly suppressed when the perfect nesting condition is broken [ 16] , 
Understanding the effect of the long-range Coulomb interaction and imperfect nesting on the 
stability of the exciton condensation and other competing phases such as a superconductivity is 
important to clarify the detailed properties of TiSe 2 and other low carrier density systems. 
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